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We present a new, universally applicable toy model of colour reconnections in hadronic final states. The 
model is based on hadronising strings and has one free parameter. We next present an implementation of 
this model in the PYTHIA event generator and provide several parameter sets ('tunes'), constrained by fits to 
Tevatron minimum-bias data. Finally, we consider the sensitivity of a simplified top mass analysis to these 
effects, in exclusive semi-leptonic top events at the Tevatron. A first attempt at isolating the genuine non- 
perturbative effects gives an estimate of order <5m top ~ ±0.5 GeV from non-perturbative uncertainties, and a 



further 5 m 



top 



±1 GeV from shower effects. 
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I. INTRODUCTION 

With increasing statistics and improved analysis tech- 
niques, a truly precise measurement of the top quark mass 
now seems feasible at the Tevatron experiments, reaching a 
final uncertainty at or below 1.5 GeV H. This is all the more 
impressive given that the top mass is a highly non-trivial ob- 
servable, involving both jets and leptons. Moreover, it fur- 
nishes an important motivation to reconsider which theoreti- 
cal aspects are relevant, at the 1 GeV level, and whether they 
are sufficiently well under control. Ultimately, this question 
will also be relevant for a range of proposed high-precision 
measurements at the LHC. 

In particular for hadronic final states, a sophisticated array 
of corrections are applied to the experimental raw data before 
the actual observable is evaluated 01 0, SI] ■ Due to the in- 
creasingly advanced procedures mandated by high precision, 
it is not straightforward to predict how uncertainties in the 
modelling affect the final answer; instead, dedicated studies 
are required to establish whether theoretical models are suf- 
ficiently well constrained and/or whether modified measure- 
ment strategies could ultimately be more fruitful. 

On the theoretical side, techniques for consistent matching 
between perturbative parton showers and fixed-order calcula- 
tions have been improved and generalised in recent years (for 
reviews see e.g. Jj, S la]), with some work focusing specif- 
ically on top production (3 HI HI OH. The structure of the 
und erly in g event (UE) has also received increasing attention 
OlESGJitll El El], with theoretical developments here fo- 
cusing on resummations of multiple perturbative interactions 
(MPI) Q3 [3 El US]. Non-perturbative aspects, on the other 
hand, still suffer from being hard to quantify, hard to test, and 
hard to calculate. In this study, we focus on one particular 
such source of uncertainty: colour reconnection effects in the 
final state. 

We begin by briefly discussing some general aspects of 
colour reconnections, including the role they already play in 
current descriptions of hadron collisions. We next present 



several explicit models, with parameters constrained by Teva- 
tron minimum-bias distributions. Finally, we apply the mod- 
els in the context of semileptonic top events at the Tevatron 
and study the sensitivity of simplified top mass estimators to 
the model variations. Some previous work leading up to this 
report can be found in I2III22I1 . 



n. COLOUR RECONNECTIONS 

In a first study of colour rearrangements, Gustafson, Pet- 
tersson, andZerwas (GPZ) 1 23] observed that, e.g. in hadronic 
WW events at LEP, colour interference effects and gluon ex- 
changes can cause 'crosstalk' between the two W systems. In 
the GPZ picture, the corresponding changes occurred already 
at the perturbative QCD level, leading to predictions of quite 
large effects. 

Sjostrand and Khoze (SK) I24I [25ll subsequently argued 
against large perturbative effects and instead considered a sce- 
nario in which reconnections occur only as part of the non- 
perturbative hadronisation pha se. Starting from the Lund 
string fragmentation model 12611 . SK argued that, if two QCD 
strings overlap in space and time, there should be a finite 
possibility for them to fuse or cut each other up (see e.g. 
B27I0 . However, since it is not known whether the QCD vac- 
uum more resembles a (chromomagnetic) Type I or Type II 
superconductor, SK presented two limiting-case models, re- 
ferred to as SK-I and SK-II, respectively. Both models re- 
sulted in effects much smaller than for GPZ, leading to a pre- 
dicted total uncertainty on the W mass from this source of 
<7 Mw < 40 MeV. SK also performed a study of QCD inter- 
connection effects in tt production [28:], but only in the con- 
text of e + e~ collisions. 

Subsequently, a number of alternative models have also 
been proposed, most notably the ones proposed by the Lund 
group, based on QCD dipoles HI S3 Hi, and one based on 
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clusters by Webber 02H . Apart from WW physics, colour re- 
connections have also been proposed to model rapidity gaps 
iSHSSSHS^l an d quarkonium production S3- 

Experimental investi gati ons at LEP did not find conclusive 
evidence of the effect ll38l SH H3, El, but were limited to 
excluding only the most dramatic scenarios, such as GPZ 
and versions of SK-I with the recoupling strength parame- 
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ter close to unity. Furthermore, in hadron collisions the ini- 
tial state contains soft colour fields with wavelengths of order 
the confinement scale. The presence of such fields, uncon- 
strained by LEP measurements, could impact in a non-trivial 
way the formation of colour strings at the time of hadroni- 
sation lf33l l34tl . And finally, the underlying event produces 
an additional amount of displaced colour charges, translating 
to a larger density of hadronising strings between the beam 
remnants. It is not known to what extent the collective hadro- 
nisation of such a system differs from a sum of independent 
string pieces. 

As the starting point for a concrete mo del we take the de- 
scription of hadron collisions developed in | 



19LM, as im- 



plemented in the Pythia event generator 14211 . In particular, 
this implies that the underlying event is obtained by a resum- 
mation of perturbative QCD 2^2 scatterings. The required 
multi-parton luminosities are obtained from the standard 1- 
parton ones, augmented with impact-parameter dependence 
and imposing flavour and momentum conservation IU9I1 . 

The system of coloured partons emerging from the pertur- 
bative phase is dual to a set of colour-singlet QCD dipoles 
ll43l l44tl . with the transformation between the two being 
uniquely determined in the large-A^ limit (modulo a small 
ambiguity between the hard scattering and underlying event 
initiators lfl9lo . In the absence of colour reconnections or 
other collective phenomena, each such dipole translates di- 
rectly to a hadronising string piece. 

Colour reconnections can then be introduced by defining a 
finite probability ^reconnect f° r eacn dipole to undergo some 
form of modification before it hadronises. We shall here de- 
fine the probability for the dipole to 'survive' as 

"int 

^kccp = 1 "^reconnect = n(!-^) = ( 1 -^)"' nt ' W 
i=l 

where £r represents an averaged probability for the 
dipole/string to interact and the scaling with the number of 
multiple interactions nj n t, is intended to give a rough rep- 
resentation of a scaling with the number of strings between 
the remnants, each of which the dipole could have interacted 
with. In principle, we could have gone further, noting that a 
dipole cannot interact with itself (n; nt — > n;„t — 1), that gluon 
exchange stretches two strings, rather than one, between the 
remnants (n; nt — > 2n g + n q ), the possibility to interact with 
the background vacuum (n- m t — > n- m t + c), etc. Given the 
uncertainties, the present scaling should be a reasonable first 
approximation, leaving the possibility open that future studies 
may require a more sophisticated behaviour. 

By consequence, for e + e~ collisions the effective recon- 
nection probability is simply while it grows from a 
minimum value close to £r for low-multiplicity (peripheral) 
hadron collisions to larger values in harder, more central col- 
lisions, such as the high-multiplicity tail of min-bias or top 
production. The model thus yields some of the expected qual- 
itative behaviour while leaving only a single free parameter, 
£r, to be determined from data. 

The dipoles which do not survive define an overall colour 
neutral system of (anti-)triplet charges (a gluon is represented 
as the sum of a triplet and an antitriplet in the N c — > 00 limit) 



for which a new string topology is to be determined. The ba- 
sis of our model is an annealing-like algorithm 112 111 which 
attempts to minimise the total potential energy, as represented 
by the string length measure A l45ll46ll . here given for mass- 
less partons for simplicity: 



N 2 
11 M? 



(2) 



where i runs over the number of string pieces, N, with invari- 
ant masses to;, and Mo is a constant of order the hadronisation 
scale. The actual measure used by the algorithm is the four- 
product of the momentum vectors of the dipole endpoints. We 
note that the minimisation of a similar measure also lies at the 
heart of an earlier model for string re-interactions proposed by 
Rathsman Il35ll . the Generalized Area Law (GAL), to which 
we plan to return in a future study. 

More aggressive models could still be constructed, e.g. by 
reducing the risk of the annealing procedure getting trapped in 
shallow local minima, but we do not consider this a critical is- 
sue. One could also be more selective about which dipoles to 
include in the annealing; here, we simply select a random set 
of 'active' dipoles, whereas in a more aggressive model one 
could have introduced a preference for dipoles which have the 
most to gain in A. Conversely, less aggressive models could be 
motivated by arguing that fairly long-lived resonances should 
be able to 'escape' the mayhem and hadronise independently. 
Presumably, this would be particularly relevant for colour sin- 
glet resonances, such as the W, and in a more sophisticated 
treatment a gradual suppression with distance from the central 
hadronising region, or, more precisely, the distance between a 
pair of interacting dipoles, would be expected. However, since 
both the top and W have decay lengths of order 0. 1 fm, for the 
present we treat their decay products as fully participating in 
the swapping of colours. 

Below, we investigate three variants of the algorithm; 
Type in which the collapse of the colour field is driven by 
free triplets only (gluons are sequentially attached to string 
systems starting from quarks), which naturally suppresses the 
formation of small closed gg string loops and is simultane- 
ously numerically the fastest, Type 1 in which gg loops are 
suppressed by brute force, and Type 2 in which gg loops are 
not suppressed. The physical question behind this issue is, 
very briefly stated, whether, at the non-perturbative transition 
level, gluons should be interpreted merely as representations 
of transverse excitations (or 'kinks') on strings whose main 
topology is defined by their quark endpoints, or whether the 
gluons should be allowed to play a more independent dynami- 
cal role. Neither is likely to be the full answer, and so by these 
variations we seek to explore some measure of the associated 
uncertainty. 

The strength parameter £r remains to be determined. In 
principle, the constraints from LEP and diffractive processes 
would be prime candidates for this task. However, since we 
are here explicitly concerned with possible breakdowns of 
jet universality, the applicability of such constraints to hard 
inelastic hadron collisions would be, at least, questionable. 
Thus, although these connections are certainly worth explor- 
ing in more depth, we note that a smaller extrapolation relative 
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to the process we are interested in can be obtained by simply 
using non-top Tevatron data. Here we consider minimum-bias 
(inelastic, non-diffractive) events at Run II of the Tevatron, the 
large-multiplicity tail of which should be fairly directly related 
to top production. In a more elaborate study, this should be 
extended to include also Drell-Yan and dijet data, particularly 
with Q 2 ~ m 2 op . 

The complex nature of hadron collisions, and in particular 
the uncertainties associated with the underlying event, imply 
that we cannot just correlate a single distribution with and 
be done with it. Instead, the colour reconnection strength must 
be determined as part of a more general fit or 'tune' to several 
minimum-bias distributions simultaneously. As a first step, 
we here use the charged hadron multiplicities, P(N c b.), and 
the mean p± as a function of multiplicity, (p±) (iV c h). 

The naive expectation from an uncorrected system of 
strings decaying to hadrons would be that (p±) should be in- 
dependent of JVch, and to first approximation equal to the LEP 
fragmentation p±_ width. (With Pythia, the best fit for the 
non-perturbative component of this is (p_i_) NP ~ 0.36 GeV, 
to give the order of magnitude.) Already at SppS, however, 
and more recently at RHIC and the Tevatron, such a constant 
behaviour has been convincingly ruled out. Currently, models 
which successfully describe the ( p±) (iVch) d istribution, such 
as R. Field's 'Tune A' and others $1% \1% 11611 . do so by incor- 
porating very strong ad hoc correlations between final-state 
partons from different interactions. We emphasise that these 
correlations are not chosen at random but are constructed to 
minimise the resulting string length, i.e. similarly to our mod- 
els here. Thus, although colour reconnections are not explic- 
itly part of these models, an implicit effect with similar conse- 
quences is still needed, at a seemingly large magnitude. This 
observation alone serves as a significant part of the motivation 
for our study. 

Given the good agreement between Tune A and Tevatron 
data, and given the difficulty in obtaining the data itself, we 
constrain the new models simply by comparing to Tune A. 
Tab. U gives a list of 8 different Pythia (v.6408) parameter 
sets, which almost unavoidably combine variations of both 
perturbative and non-perturbative aspects. However, we have 
chosen the models such that, by comparing all of them, it 
should be possible to separate the perturbative from the non- 
perturbative components, at least to a first approximation. 

Tunes A and DW both pertain to the 'old' UE model 
1 17] and are the result of careful comparisons to CDF data 
1 12, [l3i[l6ll . for Tune A including underlying event only, and 
for Tune DW also Drell-Yan data. We again emphasise the 
large role played by non-trivial colour correlations in these 
tunes, which was originally introduced to improve the fit to 
the high-p^ tail of hadron spectra J47ll . The Apt model is a 
re-tune of Tune A, with the original virtuality-ordered final- 
state showers replaced by the new p± -ordered ones l20ll (with 
Aqcd obtained from a fit to ALEPH data B). We note that, 
by default, both choices of ordering variable incorporate ma- 
trix element merging for hard jet radiation 01 for both top and 
W decays. We thus expect differences in the out-of-cone ef- 
fects from hard perturbative radiation to be small. Similarly, 
Aqr is again identical to Tune A, except that it starts from 




Tevatron min-bias - PYTHIA 6.408 
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FIG. 1: Comparison of the models/tunes discussed in the text. In- 
elastic non-diffractive (min-bias) events in pp collisions at a/s = 
1960 GeV. Top: Charged multiplicity distribution. Bottom: mean 
p± in GeV, as a function of charged multiplicity. The main point 
is not the precise predictions for each tune, but rather that they all 
roughly agree, with the notable exception of the NOCR one. 



uncorrected UE string systems and then applies the Type 
colour annealing model presented here as an afterburner. The 
S a models pertain to the 'new' interleaved UE model l2(ill . 
with p^-ordered showers and Type a colour annealing, re- 
spectively — constrained using the distributions in Fig.Q]and, 
for the FSR A QCD , ALEPH event shape data JH . The NOCR 
model is shown for reference only. It is the only one that 
does not incorporate explicit non-trivial colour correlations, 
and even then care has been taken to exploit the initial-state 
colour ambiguity mentioned above to the fullest. 

In Fig. Q] the 8 models in Tab. U are compared on the 
iV c h and (p±) (N^) distributions. A good description of the 



4 



Parameter (PYTHIA V.6408+) 


DW A Apt Aqr 


So Si S 2 NOCR 


UE model MSTP(81) 
UE infrared regularisation scale (at y/s = 1800 GeV) PARP(82) 

-"- Qpnlino nnwpr with . /q PARPiQfli 

UE hadron transverse mass distribution MSTP(82) 
-"- parameter 1 PARP(83) 
-"- parameter 2 PARP(84) 

UE total gg fraction PARP(86) 


1 ('old' [IT]) 
1.9 2.0 2.1 2.0 

25 ffa<;r"> 
4 ('double Gaussian') 
0.5 
0.4 

1.0 0.95 0.95 0.66 


21 ('new' [20]) 
1.85 2.1 1.9 2.05 
1 6 r'slnw"! 
5 ('ExpOfPow') 
1.6 1.4 1.2 1.8 
n/a 
n/a 


ISR infrared cutoff PARP(Wl 

lulV 11111 ill ^LIlVJll I ./\.IVI \ } 

ISR renormalisation scale prefactor PARP(64) 
ISR Q^ax factor PARP(67) 
ISR infrared regularisation scheme MSTP(70) 
ISR FSR off ISR scheme MSTP(72) 


1 25 1 10 10 
0.2 1.0 1.0 1.0 
2.5 4.0 4.0 4.0 

n/a 
n/a 


C= PARPf821 1 
1.0 
n/a 

2 2 2 
10 


FSR model MSTJ(41) 
FSR Aqc D PARJ(81) 


2 2 12 2 
0.29 0.29 0.14 0.29 


(p± -ordered) 
0.14 


BR colour scheme MSTP(89) 
BR composite x enhancement factor PARP(79) 
BR primordial k T width (| k T | > PARP(9 1 ) 
BR primordial k T UV cutoff PARP(93) 


n/a 
n/a 

2.1 1.0 1.0 1.0 
15.0 5.0 5.0 5.0 


1112 

2 2 2 3 
n/a 
5.0 


CR model MSTP(95) 
CR strength § B PARP(78) 
CR gg fraction (old model) PARP(85) 


n/a 6 
n/a 0.25 
1.0 0.9 0.9 0.0 


6 2 4 1 
0.2 0.35 0.15 0.0 
n/a 



TABLE I: PYTHIA parameters |43l . divided into a few main categories: UE (underlying event), ISR (initial state radiation), FSR (final-state 
radiation), BR (beam remnants), and CR (colour reconnections). The UE reference energy for all models is parp (8 9) =1800 GeV, and all 
dimensionful parameters are given in units of GeV. MSTP ( 95 ) =2 , 4 , 6 corresponds to CR types 1, 2, and 0, respectively, in the text. 



charged multiplicity distribution is obtained in all cases. Si- 
multaneous good agreement with (p±) (N^) is only obtained 
for the models incorporating non-trivial colour correlations - 
the notable exception being the NOCR model which exhibits 
the close-to-constant behaviour of uncorrelated string decays 
discussed above. We interpret this behaviour as represent- 
ing a concrete example of a data-driven motivation to develop 
ideas of hadronisation beyond the current cluster/string mod- 
els, which have remained essentially frozen since the LEP era. 
This is not to imply that these models are intrinsically unsat- 
isfactory or that hadronisation at LEP should be perceived as 
being a completely separate story, but only that 1) we expect 
that the many recent improvements on the perturbative side 
imply that there is less 'wriggle' room for the non-perturbative 
physics, and hence the latter could presumably be better con- 
strained today than a decade ago, and 2) the size of possible 
differences between non-perturbative effects in different envi- 
ronments should be more fully explored, where current mod- 
els are normally limited to the assumption of jet universality. 

Returning to the question at hand, the remaining models in 
Tab. U all describe the two distributions in Fig. Q] within an 
acceptable margin, at least as gauged by the spread between 
the two more elaborate Tevatron tunes, A and DW. 



III. THE TOP MASS AT THE TEVATRON 

Assuming the in situ extrapolation discussed above to be at 
least moderately reliable, we now apply the same models and 
parameters in the context of top production at the Tevatron. 



More specifically, we concentrate on semileptonic ti events, 
i.e. ti — > bbqqiv. As a first step towards estimating the sensi- 
tivity of experimental top mass observables, we consider the 
impact on a simplified measurement which roughly approxi- 
mates the key ingredients used in current Tevatron analyses. 
As already mentioned, it is not possible to cleanly separate 
the CR effect from other sources of variation in these models. 
Nonetheless, by grouping models with similar parton show- 
ers and studying both the variations within and between the 
groups, we are able to make some headway. Differences be- 
tween the results obtained using the same parton shower but 
different CR models for the same generated top mass may then 
be interpreted as a first estimate of the uncertainty on mt op 
due to genuine non-perturbative effects, while the uncertainty 
between groups with different shower models is interpreted as 
having a perturbative origin. 



A. Real Top Mass Measurement 

To set the stage for the subsequent analysis, let us first sum- 
marise the methods used in actual top mass reconstructions at 
the Tevatron, focusing on the ingredients they employ to go 
from detector level events to a reconstructed top mass. The 
relevant measurements performed in the semileptonic channel 
by the CDF and D0 collaborations 61 |U HI HI HI use 
one of three methods. 

The template method compares the distribution of kinemat- 
ically reconstructed top mass values to templates obtained for 
various nominal top mass values from simulation with full 
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detector description (including background). The matrix ele- 
ment method computes the event-by-event likelihood that the 
observed kinematic configurations stem from events of a given 
top mass. Maximising the total likelihood of the observed 
sample yields the final result. Finally, the ideogram method 
reconstructs top mass values in each event and then builds the 
likelihood of observing that value with the given resolution 
as a function of the true value. Again, maximising the total 
likelihood of the observed sample yields the final result. 

All three methods are based on so-called reconstructed 
physics objects, i.e. jets, identified charged leptons, and miss- 
ing transverse energy, which fulfil certain selection criteria. 
For the methods with explicit reconstruction of the top mass 
an assignment of jets to one of the two top quarks has to be 
made. Often constrained fits are used to improve the experi- 
mental resolution by requiring that the two reconstructed top 
masses are equal and that the known VK-mass is reconstructed 
in the decay of each of the top quarks. In this case an assign- 
ment of each jet to the VK-boson or 6-quark within each top 
decay product is also required. 

The methods are calibrated, i.e. the performance of a given 
implementation is measured on fully simulated events. Any 
deviation of the reconstructed from the generated top masses 
is corrected for; this implies that the top mass measurement is 
limited by the precision of the simulation used for the calibra- 
tion. 

Since summer 2005 all methods have been extended to 
tackle the dominant experimental systematic uncertainty, the 
jet energy scale (JES), by simultaneously fitting the top mass 
and the jet energy scale, with the additional constraint that the 
reconstructed mass of the hadronically decaying W present in 
each event should be consistent with the known mw- 

The main features of a top mass measurement are thus util- 
isation of reconstructed physics objects, assignment of each 
jet to a specific top decay product, correction for an overall 
JES factor, and calibration of the method. We shall now con- 
struct a simplified top mass estimator which embodies these 
four ingredients. 



B. Toy Top Mass Measurement for Generator Level 

To obtain an estimate of the influence of the above CR/UE 
models on top mass measurements a toy mass measurement 
on events generated by PYTHIA was implemented. The study 
is performed at the generator level, here meaning after hadro- 
nisation and hadron decays but without detector simulation. 

In this simplified analysis electrons and neutrinos are 'iden- 
tified' by looking at the generator truth. Jets are reconstructed 
on final-state particles using cone jets of AR = 0.5 fsil l55ll 
and are required to have apr > 15GeV. Only semileptonic 
events with exactly four reconstructed jets are considered for 
further analysis. 

Jets are assigned to the top decay parton with the lowest 
AR distance. Events without a unique one-to-one assignment 
are discarded. In the samples used here, between 25% and 
27% of the semileptonic events fulfil these requirements. 

The top mass is now computed event by event from the sum 



of the four-momenta of the three jets from the hadronically 
decaying top in the event. The top mass for the full sample is 
then obtained by fitting a Gaussian to the distribution of recon- 
structed masses. To minimise the importance of the tails of the 
distribution, the fit range is restricted to a window of ±15 GeV 
around the top mass. This width corresponds approximately 
to the experimental resolution observed in current measure- 
ments. The fit is iterated until it settles on a range symmetric 
around the final result, m^ p . 

The same method can also be applied for measuring the 
VF-mass from the two jets assigned to the partons from 
the W decay. The resulting VF-mass measurement can 
be used to compute a JES correction factor: sjes = 
80.4GeV/mvK- The fitted top mass above can then be cor- 
rected by this factor, to produce a JES-scaled top mass value, 

^scaled _ imfit 

m to P = S JES • m top . 

This full procedure is repeated for several different values 
of the generated top mass, between 165 and 185 GeV, in steps 
of 1 GeV. Both the scaled and unsealed fitted top masses ex- 
hibit a completely linear dependence on the input top mass. 
Observed slopes are consistent with 1 at the 2% level. This 
indicates that the fit procedure is indeed stable and has the 
desired dependence on the physical quantity, cf. Fig. [2] 



C. Differences from different Models 

In real mass measurements, the offset and slope of the 
straight line that describes the reconstructed vs. the generated 
top mass is used to calibrate the given top mass procedure. 
However, this calibration must necessarily use only one spe- 
cific CR/UE model. By virtue of the tuning we performed, 
there is a genuine ambiguity of which model to choose. Dif- 
ferences between the individual model calibrations therefore 
lead to uncertainties on the top mass results. As an example, 
Fig.|2]shows the calibration curve obtained for Tune A before 
JES rescaling. 

Fig. |3] summarises our central results. It shows the offsets 
before (left) and after (right) scaling with the JES correction 
factor. The offsets, Amtop, for each model are obtained from 
a straight line fit to the calibration curve evaluated at mf° p = 
175 GeV, and have statistical precisions (determined from the 
spread of the data points) of ~ 0.1 GeV. 

The top masses in the uncorrected fits (dots, left column) 
come out somewhat lower than the input mass, principally due 
to out-of-cone corrections. Including the JES correction, i.e. 
scaling all jets by the factor necessary to get the right hadronic 
W mass, the points move to the right (squares), even to the 
point of over-correcting the top mass. 

Again, our central point is that, while for any particular 
model a further, constant offset would be sufficient to cali- 
brate the measurement to coincide with the input mass, the 
spread between models cannot be dealt with in this way. It is 
the ambiguity coming from not knowing which offset value to 
correct for that we interpret as the uncertainty on the top mass. 

It therefore seems significant that the various models ex- 
hibit differences of about ± 1 . 1 GeV and ±1.5 GeV for the off- 
sets of m^op and m^p lcd , respectively. Explicit checks vary- 
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FIG. 2: Calibration curve obtained for Tune A, before JES rescaling. 
A similar plot was made for each model in Tab. U and their rela- 
tive offsets compared, both before and after JES rescaling. The inset 
shows the Gaussian fit to the distribution reconstructed top masses 
from the hadronic event side for the specific point mf™ = 175 GeV. 



-10 




Tune A 
Tune A-CR 
Tune A-PT 
Tune DW 
Tune BW 

SO 
SI 
S2 

NoCR 

I i i 

-5 



Am, 



FIG. 3: Comparison of calibration offsets obtained for each model, 
in GeV, here including an additional parameter set, 'BW, from Rick 
Field. On the left are the results obtained before JES rescaling (dots) 
and on the right after rescaling (squares). The coloured bands group 
models with the same final-state shower (green: virtuality-ordered, 
blue: p± -ordered). The statistical precision due to the finite number 
of generated events is at the 0.1 GeV level. 



ing both the fit range and fit function produced variations no 
larger than ~ 20% in these numbers, hence at this level the 
effect appears genuine. Without additional constraints from 
data, it translates directly into an uncertainty on the recon- 
structed top mass. 

To extricate the genuinely non-perturbative part of this, we 
note that the models fall into two broad classes: those that 
utilise the "old" virtuality-ordered final-state parton shower 
and those that utilise the "new" p± -ordered one. The largest 
component of the difference is between these two classes, 
hinting at a perturbative origin for most of it, which, at least 
to some extent, should already be present in Tevatron analyses 
via the Pythia-Herwig systematic. 

Within each class, we still observe differences roughly of 
order ±0.5 GeV on the top mass, which we are more confi- 
dent in assigning a non-perturbative origin. Note, however, 
that this still lumps genuine CR effects together with other in- 
frared ambiguities, such as infrared regularisation and renor- 
malisation procedures, the treatment of beam remnants, etc. 

Real mass analyses may have a different sensitivity to the 
model differences. The size of the effect in this first study, 
however, suggest a need for further in depth analyses. If the 
sensitivity we observe here is confirmed for real mass anal- 
yses, we hope the question may be turned around, and that 
in situ measurements can be used to gain further information 
about the interesting physics effects that may be present. 



IV. SUMMARY AND CONCLUSION 

We have presented a set of new, universally applicable mod- 
els to study colour reconnection (CR) effects in hadronic final 
states. The models are based on hadronising strings and ap- 
ply an annealing-like algorithm to minimise a measure of the 
classical potential energy, with a freely varying strength pa- 
rameter running from zero to unity. A scaling is included such 
that the survival probability of a given string piece decreases 
as a function of the number of perturbative scatterings in the 
underlying event. The models are implemented in a publicly 
available version (v.6408+) of the PYTHIA event generator. 

To constrain the CR strength parameter we have used 
Tevatron minimum-bias distributions, specifically P(iV c h) 
and (p±) (iVch). Taking the results obtained with the CDF 
'Tune A' as a benchmark, we present several alternative pa- 
rameter sets exploring the possible combinations of show- 
ers, underlying-event modelling, and colour reconnections. 
As a further, data-driven motivation, we argue that current 
underlying-event descriptions, including 'Tune A', already in- 
clude strong non-trivial colour correlations. 

As a first application, we have investigated the influence 
of changing the underlying physics model, including CR, UE, 
and shower effects, on a range of simplified Tevatron top mass 
measurements. The models we consider exhibit individual 
variations of about ±1.5 GeV on the reconstructed top mass. 
While this is comparable to systematic uncertainties quoted 
for present top mass measurements, it has so far only partly 
been considered in the current analyses. 

Of the total variation we attribute about ±lGeV to per- 
turbative effects and about ±0.5 GeV to non-perturbative 
sources. 
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Our conclusion for the present is thus twofold: firstly, 
colour reconnections in hadron collisions appear to be a both 
experimentally and theoretically motivated possibility, one 
which should be explored as part of developing a more de- 
tailed picture of hadron collisions. Secondly, it appears that 
non-perturbative uncertainties, among which colour recon- 
nections hold a prominent place, are likely to be relevant in 
the drive towards sub-GeV uncertainties on the top mass at 
the Tevatron. 

It is important to now verify the size of the observed uncer- 
tainties in real mass measurements and to increase the amount 
of non-top data used to constrain the models; Drell-Yan pro- 
duction, and in particular its high-mass tail, is likely to be 
useful in reducing the initial-state shower ambiguities, while 
the infrared effects could be further probed by expanding on 
the number of minimum-bias distributions, as well as includ- 
ing underlying-event studies in dijet and Drell-Yan produc- 
tion, again in particular when Q ~ m top where the required 
extrapolation is presumably minimal. The connection with 
the LEP data (and, possibly, diffractive physics) should also 
be explored, although as we have noted the assumption of jet 



universality should probably not be treated as inviolate in this 
context. 

On the theoretical side, we hope that the arguments we have 
presented will stimulate curiosity, and eventually activity, in 
this now somewhat dormant field. Along the intersection of 
the two communities, it would be interesting to explore alter- 
native measurement strategies and, as a last resort, a combined 
tuning and top mass fit. A final follow-up we envision is to ex- 
trapolate in energy to evaluate the impact on precision studies 
at the LHC. 



Acknowledgements: 

We are grateful to F. Canelli, F. Fiedler, K. Hamacher, 
J. Rathsman, M. Sandhoff, and T. Sjostrand for enlightening 
discussions and valuable comments on the draft. P.S. is sup- 
ported by Fermi Research Alliance, LLC, under Contract No. 
DE-AC02-07CH1 1359 with the United States Department of 
Energy 



[1] Tevatron Electroweak Working Group (2006), 
http : / /www-cdf . f nal . gov/physics/new/top/top .html, 
see also hep-ex/0608032. 

[2] A. Bhatti et al„ Nucl. Instrum. Meth. A 566, 375 (2006), hep- 
ex/05 10047. 

[3] J. D'Hondt, S. Lowette, J. Heyninck, and S. Kasselmann 

(2006), CERN-CMS-NOTE-2006-025. 
[4] M. L. Mangano and T. J. Stelzer, Ann. Rev. Nucl. Part. Sci. 55, 

555 (2005). 
[5] S. Hoche et al. (2006), hep-ph/060203 1 . 
[6] T. Sjostrand (2006), hep-ph/061 1247. 

[7] E. Norrbin and T. Sjostrand, Nucl. Phys. B 603, 297 (2001), 

hep-ph/0010012. 
[8] S. Frixione, P. Nason, and B. R. Webber, JHEP 08, 007 (2003), 

hep-ph/0305252. 
[9] K. Hamilton and P. Richardson (2006), hep-ph/0612236. 
[10] M. L. Mangano, M. Moretti, F. Piccinini, and M. Treccani, 

JHEP 01, 013 (2007), hep-ph/061 1 129. 
[11] T. Affolder et al. (CDF), Phys. Rev. D 65, 092002 (2002). 
[12] R. Field and R. C. Group (2005), hep-ph/05 10198. 
[13] R. Field (CDF), AIP Conf. Proc. 828, 163 (2006). 
[14] S. Abdullin et al. (2006), CERN-CMS-NOTE-2006-033. 
[15] D. Acosta et al. (2006), CERN-CMS-NOTE-2006-067. 
[16] M. Albrow et al. (2006), hep-ph/0610012. 
[17] T. Sjostrand and M. van Zijl, Phys. Rev. D 36, 2019 (1987). 
[18] J. M. Butterworth, J. R. Forshaw, and M. H. Seymour, Z. Phys. 

C 72, 637 (1996), hep-ph/9601371. 
[19] T. Sjostrand and P. Z. Skands, JHEP 03, 053 (2004), hep- 

ph/0402078. 

[20] T. Sjostrand and P. Z. Skands, Eur. Phys. J. C 39, 129 (2005), 
hep-ph/0408302. 

[21] M. Sandhoff and P. Skands (2005), FERMILAB-CONF-05- 
518-T. In Les Houches 'Physics at TeV Colliders' 2005 SM 
and Higgs Working Group: Summary report,/hep-ph 0602198. 

[22] M. Sandhoff (2005), Diploma Thesis, Bergische Universitat 
Wuppertal, WU D 05-08. 

[23] G. Gustafson, U. Pettersson, and P. M. Zerwas, Phys. Lett. B 



209, 90(1988). 

[24] T. Sjostrand and V. A. Khoze, Phys. Rev. Lett. 72, 28 (1994), 
hep-ph/93 10276. 

[25] T. Sjostrand and V. A. Khoze, Z. Phys. C 62, 281 (1994), hep- 
ph/93 10242. 

[26] B. Andersson, G. Gustafson, G. Ingelman, and T. Sjostrand, 

Phys. Rept. 97, 31 (1983). 
[27] X. Artru, Phys. Rept. 97, 147 (1983). 

[28] V. A. Khoze and T. Sjostrand, Phys. Lett. B 328, 466 (1994), 

hep-ph/9403394. 
[29] G. Gustafson and J. Hakkinen, Z. Phys. C 64, 659 (1994). 
[30] L. Lonnblad, Z. Phys. C 70, 107 (1996). 

[31] C. Friberg, G. Gustafson, and J. Hakkinen, Nucl. Phys. B 490, 

289 (1997), hep-ph/9604347. 
[32] B. R. Webber, J. Phys. G 24, 287 (1998), hep-ph/9708463. 
[33] W. Buchmuller and A. Hebecker, Phys. Lett. B 355, 573 (1995), 

hep-ph/9504374. 

[34] A. Edin, G. Ingelman, and J. Rathsman, Phys. Lett. B 366, 371 

(1996) , hep-ph/9508386. 

[35] J. Rathsman, Phys. Lett. B 452, 364 (1999), hep-ph/9812423. 
[36] R. Enberg, G. Ingelman, and N. Timneanu, Phys. Rev. D 64, 

114015 (2001), hep-ph/0 106246. 
[37] A. Edin, G. Ingelman, and J. Rathsman, Phys. Rev. D 56, 7317 

(1997) , hep-ph/9705311. 

[38] G. Abbiendi et al. (OPAL), Phys. Lett. B 453, 153 (1999), hep- 
ex/9901019. 

[39] G. Abbiendi et al. (OPAL) (2005), hep-ex/0508062. 
[40] S. Schael et al. (ALEPH), Eur. Phys. J. C 48, 685 (2006), hep- 
ex/0604042. 

[41] J. Abdallah et al. (DELPHI), Phys. Lett. B 643, 147 (2006), 
hep-ex/06 10031. 

[42] T. Sjostrand, S. Mrenna, and P. Skands, JHEP 05, 026 
(2006), http :/ /pro jects . hepf orge . org/pythia6/, hep- 
ph/0603175. 

[43] G. Gustafson, Phys. Lett. B 175, 453 (1986). 

[44] G. Gustafson and U. Pettersson, Nucl. Phys. B 306, 746 (1988). 

[45] B. Andersson, G. Gustafson, and B. Soderberg, Z. Phys. C 20, 



8 



317 (1983). 

[46] T. Sjostrand and P. Z. Skands, Nucl. Phys. B 659, 243 (2003), 

hep-ph/0212264. 
[47] R. D. Field (2007), private communication. 
[48] G. Rudolph, private communication. 

[49] A. Abulencia et al. (CDF), Phys. Rev. D 73, 032003 (2006), 
hep-ex/05 10048. 

[50] A. Abulencia et al. (CDF), Phys. Rev. D 73, 092002 (2006), 
hep-ex/05 12009. 



[51] V. M. Abazov (DO), Phys. Rev. D 74, 092005 (2006), hep- 

ex/0609053. 
[52] DO (2005), dOnote 4728.conf. 

[53] V. Abazov (DO) (2007), submitted to PRD, hep-ex/0702018. 

[54] J. M. Butterworth, J. P. Couchman, B. E. Cox, and B. M. Waugh 
(2002), MAN/HEP/2002/02, UCL/HEP 2002-02. 

[55] W. G. Piano, University of Manchester. Private communica- 
tions. 



